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Abstract 

We investigate the spectral properties of a generalized GOE (Gaussian Orthogonal En- 
semble) capable of describing critical statistics. The joint distribution of eigenvalues of 
this model is expressed as the diagonal element of the density matrix of a gas of particles 
governed by the Calogero-Sutherland Hamiltonian (C-S). Taking advantage of the corre- 
spondence between C-S particles and eigenvalues, we show that the number variance of 
our random matrix model is asymptotically linear with a slope depending on the param- 
eters of the model. Such linear behavior is a signature of critical statistics. This random 
matrix model may be relevant for the description of spectral correlations of complex quan- 
tum systems with a self-similar/fractal Poincare section of its classical counterpart. This 
is shown in detail for two examples: the anisotropic Kepler problem and a kicked particle 
in a well potential. In both cases the number variance and the A3-statistic is accurately 
described by our analytical results. 

PACS: 11.30.Rd, 12.39.Fe, 12.38.Lg, 71.30.+h 
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1 Introduction 



Random matrix ensembles (RME) are an invaluable tool in describing the level statistics 
of complex quantum systems. Typically, their range of applicability for disordered systems 
is determined by the Thouless energy which, in the metallic phase, is much larger than the 
average level spacing. In the neighborhood of a localization-delocalization transition the 
Thouless energy is of the order of the average level spacing and the wavefunctions become 
multifractal. The usual random matrix ensembles are no longer applicable. Recently, new 
random matrix ensembles [1, 2, 3, 4, 5, 6, 7, 8] depending on additional parameters have 
been proposed to describe spectral correlations in this critical case. These new models for 
critical statistics have been successfully utilized to describe the spectral correlations of a 
disordered system at the Anderson transition in three dimensions [9, 10], two dimensional 
Dirac fermions in a random potential [11], the quantum Hall transition [12] and of the 
QCD Dirac operator in a liquid of instantons [13, 7]. 

There are two different types of models for critical statistics. In the first one, devi- 
ations from Wigner-Dyson statistics are obtained by adding a symmetry breaking term 
to the GUE [2, 7]. The model is solved by mapping it to a non-interacting Fermi gas 
of eigenvalues. The second one [4] makes use of soft confining potentials and is solved 
exactly by means of g-orthogonal polynomials. Both models lead to the same spectral 
kernel for small deviations from the GUE. Based on this observation it was conjectured 
[1] that critical statistics is universal. However, the origin of the critical kernel is dif- 
ferent in both cases. In models based on a soft confining potential the critical kernel 
is obtained from a nontrivial unfolding. In models with an explicit symmetry breaking 
term, deviations from Wigner-Dyson statistics arise because the long range correlations 
between the eigenvalues are exponentially suppressed [14]. We remark that the extent of 
universality in critical statistics is still under debate. For instance, for chiral ensembles, 
the correlation functions of a model based on a symmetry breaking term [7] and a model 
with a soft confining potential [15] are different. Only the former one reproduces critical 
statistics. 

Critical random matrix models for orthogonal and symplectic ensembles have recently 
[16] been reported in the literature. Tsvelik and Kravtsov have obtained asymptotic 
expressions for the critical two level spectral function from a generalized ensemble of 
random banded matrices. An exact expression for the two level critical spectral function 
was conjectured in [16] for orthogonal and symplectic ensembles. In the context of the 
Anderson model, a similar result was conjectured by Nishigaki [10]. 

In order to describe the spectral correlations of certain pseudo-integrable billiards with 
dynamics intermediate between chaotic and integrable, Bogomolny and coworkers [17] 
have introduced a short range plasma model that interpolates between Poisson statistics 
and Wigner-Dyson statistics. The joint distribution of eigenvalues in [17] is given by the 
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classical Dyson gas with the logarithmic pairwise interaction restricted to a finite number 
k of nearest neighbors. Analytical solutions are available for general k and symmetry 
class. It turns out that this short-range plasma model reproduces the typical charac- 
teristics of critical statistics like a linear number variance, with a slope depending on k 
and, asymptotically, an exponential decay of the nearest neighbor spacing distribution. 
However the two models models are not identical. In critical random matrix models based 
on a symmetry breaking term, the joint distribution of eigenvalue can be considered as 
an ensemble of free particles at finite temperature with a nontrivial statistical interac- 
tion. The statistical interaction resembles the Vandermonde determinant, and the effect 
of finite temperature is to suppress the correlations of distant eigenvalues. In [17] this 
suppression is abrupt, in contrast to critical statistics, where the effect of the temperature 
is smooth. For further details we refer to [14]. 

In this paper we introduce a generalized GOE based on the addition of a symmetry 
breaking term to an invariant Gaussian probability distribution along the lines of [2] for 
the GUE. We show that the joint eigenvalue distribution of this model coincides with the 
diagonal element of the density matrix of a gas of particles governed by the Calogero- 
Sutherland (C-S) Hamiltonian. Using this identification we calculate the asymptotic 
behavior of the number variance from the susceptibility of the C-S partition function. 
Because the Itzykson-Zuber integral for /3 — 1 is unknown, a direct calculation of the 
correlation functions is not possible. To obtain analytical results we invoke the Kravtsov- 
Tsvelik conjecture, which states that the finite temperature modifications of the corre- 
lation functions arise only through the known finite temperature modifications of the 
kernel for (3 = 2. The validity of this conjecture is tested in two different ways. First, 
we show that it is in agreement with a conformal calculation of the asymptotic behavior 
of the two point correlation function. Second, the asymptotic behavior of the number 
variance according to the Kravtsov-Tsvelik conjecture agrees with the behavior of the 
susceptibility in the grand canonical ensemble. One of the main aims of this article is to 
show that critical statistics describes the spectral correlations of time-reversal invariant 
quantum systems with a corresponding classical phase space that has a global self-similar, 
fractal structure. This is shown in two examples, a kicked particle in a potential well and 
the anisotropic Kepler problem, by comparing the two-point level correlations with the 
analytical formula of our critical Random Matrix Model. 

The critical Random Matrix Model and its relations with the C-S model is discussed in 
section 2. In section 3, we review the Kravtsov-Tsvelik conjecture for the density-density 
correlation function of the C-S model in the low temperature limit. The validity of this 
conjecture is discussed in section 4. In section 5 we show that the level correlations of a 
kicked particle in a potential well and of the anisotropic Kepler problem are described by 
the Kravtsov-Tsvelik conjecture. Concluding remarks are made in section 6. 
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2 Definition of the Model 



A Random Matrix Model for Hermitian matrices with critical eigenvalue statistics was 
introduced in [2]. Although it is straightforward to generalize this model to the class 
of the Gaussian Orthogonal Ensemble, the absence of an explicit result for the integral 
over orthogonal matrices makes its analysis far more complicated. The model we study 
is defined by the joint probability distribution 

P(S, b)= J dMe ~^ ssT e'^ M ^ M ^ T . (1) 

Here, the N x N matrices S and M are real symmetric and orthogonal, respectively, and 
the integration measure dM is the Haar measure. From the invariance of dM it follows 
that P(S, b) is a function of the eigenvalues of S only. If the eigenvalues of S are denoted 
by x k the joint eigenvalue distribution is given by 

p(x ir --,x N ) = A({x k }) J dMe-^ 2b+1) ^ xl+b ^ M ^ XkXl - (2) 

where the Vandermonde determinant is defined by 

A({x k }) = H(x k -xi). (3) 

k<i 

Let us now consider the harmonic oscillator Hamiltonian 

H = -V| + ^cu 2 TrS 2 (4) 
where the Laplacian for symmetric matrices S is given by 

2 N d 2 1 N d 2 
i= i ab u 1 i<j ab ij 

This Hamiltonian is the sum of N(N + l)/2 independent harmonic oscillators with imag- 
inary time propagator given by [18] 

/ , , \ N(N+l)/4 

(S\e- TH \S')= - ] e"4ifc[(^ 2 +™ ' 2 )cosho;r-2TVS5'] /gx 

\47rsinhcur/ 

Since the Laplacian, V|, is the sum of a radial piece, depending only on the eigenvalues of 
S, and an angular piece, depending only on the orthogonal matrix M s that diagonalizes 

s, 
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the matrix element (S\e~ TH \S') factorizes into a radial piece and an angular piece. After 
integration over the angular degrees of freedom and putting the eigenvalues of S' equal 
to the eigenvalues of S we obtain 

(x 1 ,---,x N \A l ^({x k })e- TH ^A- 1 / 2 ({x k })\x 1 ,---,x N ) =C j dMe'^r^ ! -w-TiSmsm^ 

(8) 

where the integral over the angular matrix element has been absorbed in the normalization 
constant C. If we make the identification 

w , wcoshcjT 

= 2b and — — = 26 + 1, (9) 



sinh cut sinh cut 

the r.h.s. of this equation is exactly the joint probability distribution (2). We thus have 
shown that the joint probability distribution of our model is given by the diagonal matrix 
element of the density matrix of the Hamiltonian 

H = A^ 2 {{x k })H tad A-^ 2 ({x k }). (10) 

Using the identity, 

A-(M> g ^A(M)^(M) ^ ^ + - g J — w , (11) 
we find the Hamiltonian 

H = ~^~dx]~~A^.{x t -x j y + T^ 1 (12) 

This Hamiltonian corresponds to the Calogero-Sutherland model [19, 20] 

x — -v d 2 A . A . x — > 1 co 2 x — -v o 

»« = -&+ 5 ( r i)fa+ T K. (13) 

3 3 f=F3 v ■> ' 3 

with A = 1 and fermionic boundary conditions. We have thus shown that the joint 
eigenvalue distribution of the model (1) is given by the diagonal matrix elements of the 
./V-particle density matrix of the Calogero-Sutherland model at an inverse temperature r 
given by (9). 

The normalized eigenfunctions of the Calogero-Sutherland Hamiltonian (12) can be 
labeled in terms of the partitions of integers denoted by k (see next section). If A K 
and ty K (xi, ■ ■ ■ , xn, oj) are the eigenvalues and eigenfunctions of the C-S Hamiltonian, 
respectively, the joint eigenvalue probability distribution is given by 

p(x x , ...,x N ) = C'J2 e~ x * T ^ K (x, lu)^ k (x, to) (14) 
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where C is a constant and x — x±, • • • , x n . The ^ K (x, u) can be expressed in terms of the 
generalized Hermite polynomials [30] 

* k (x,u) = -^= e -i w S**2A 1 / 2 ({x,})ff Jfc (x v /^72,2) (15) 

where N K is a normalization constant and the eigenvalue is given by 

\ K = cu\k\. (16) 

We point out that the above relation between the Euclidean propagator in symmetric 
spaces and the C-S model at finite temperature can be extended to all nine other sym- 
metry classes in the Cartan classification of large families of symmetric spaces [21]. In 
essence, the radial part of the Laplacian in the symmetric space corresponds to a C-S 
type Hamiltonian. For a classification of C-S Hamiltonians based on the symmetry class 
we refer to [22]. 

Finally, let us mention that the interpolating role (between RMT and the Poisson 
ensemble) of b can be inferred directly from (1). Using the invariance of the measure, the 
integral over M can be replaced by an integral over the eigenvalues of M. For b — > oo, 
this partition function is dominated by matrices S that commute with arbitrary diagonal 
orthogonal matrices. This set of matrices is the ensemble of diagonal symmetric matrices 
also known as the Poisson ensemble with uncorrected eigenvalues. Critical statistics is 
obtained in the thermodynamic limit if the parameter b is scaled as 

b = h 2 N 2 . (17) 

Wigner-Dyson statistics is found for a weaker iV-dependence of b whereas a stronger N- 
dependence leads to Poisson statistics. This transition can also be understood in terms 
of the C-S model at finite temperature. At zero temperature, the probability density of 
the ground state of the C-S model (12) coincides with the joint probability distribution of 
the Gaussian Orthogonal Ensemble. In the high temperature limit, r — > 0, the positions 
of the particles become uncorrelated and the statistics of the associated matrix model is 
Poisson. 

To recapitulate, we have traded the problem of performing an integral over the or- 
thogonal group by the physical task of finding the diagonal element of the density matrix 
of an ensemble of particles governed by the C-S Hamiltonian. 

2.1 Excited states and Zonal Polynomials 

In this section we discuss explicit solutions of the excited eigenfunctions of the C-S Hamil- 
tonian (12) and argue to what extent they are useful for the evaluation of correlation 
functions from the joint eigenvalue distribution (2). 
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The probability density of the ground state of the C-S model (12) coincides with the 
joint probability distribution of the Gaussian Orthogonal Ensemble [19]. This observation 
together with the conjecture of the solvability of the C-S model was already made in the 
pioneering articles of Calogero [20] and Sutherland [19]. Later, Sutherland [23] obtained 
a non-orthogonal set of solutions. The problem of finding a set of orthogonal solutions 
for these excited states was recently solved by Forrester, Ha and Serban [24, 25, 26] 
who expressed the wavefunctions of the excited states in terms of the symmetric Jack 
polynomials [27]. For the special values of the coupling constant related to GOE and 
GSE the Jack polynomials have a geometrical interpretation and are usually called zonal 
polynomials [28]. Unfortunately, there is no closed formula neither for the Jack nor 
for the zonal polynomials [28]. Since explicit calculations rely on recurrence relations, 
numerical work is needed to evaluate polynomials of high degree. In our case, due to the 
harmonic potential, the excites states are given by the generalized Hermite (or Hidden- 
Jack) polynomials [29] which can be expressed in terms of Jack polynomials [29, 30]. 

The generalized Hermite polynomials in (14) can be expressed in terms of zonal poly- 
nomials by means of a Mehler type formula [30]. 

p(x u ...,x N )cx. e^^ x lA({x k }) ]T ^H K (x^J/2, 2)H K (x^J/2, 2) 



oc A({x k })e-^ coth ( T ^ x * V — " V — -(18) 



where the C^(x) are the symmetric Jack polynomials as defined in [30], x — xi, . . . , x n 
and k labels the partitions of the integers (there is a polynomial for each partition) and 
the sum runs over all the partitions. Furthermore, r and u> are related to b through 
equation (9). 
The kernel 

F ° (x ' y) = Vw cP(i«) (19) 

has been studied extensively [28, 30]. Using the result for the r — > limit of the kernel 
[30], 

1 N 

F (x/v^, y/Vr) oc J] c™ iw/2 (20) 



one easily shows that the eigenvalues of our model are uncorrelated in the high tempera- 
ture limit. 

In the zero temperature limit (GOE), the joint distribution can be represented as a 
quaternionic determinant and the integrations can be performed by means of a "kernel 
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relation". At nonzero temperatures, the kernel F (x, y) satisfies the "kernel relation" [30], 
J dn(y)F (2y,z)F (2y,x) oc F (2z, x )e~^T,M+^) j (21) 

where dp(y) = TiiLi e~% y i A({y k })dyi . . . dy^, but it is not known 1 whether the joint 
eigenvalue distribution can be expressed in terms of quaternionic determinants. Further 
progress could rely on exploiting the orthogonality relations [28] verified by the Zonal 
polynomials. If proceeding so, the partition function associated with (1) can be evaluated 
exactly. As expected, it coincides with the one encountered for particles obeying fractional 
statistics in one dimension [31]. Because of technical problems, the same strategy fails for 
the spectral density and higher order spectral correlation functions. 

For the special case of iV = 2, an explicit calculation of the joint distribution of 
eigenvalues (2) coincides with the result obtained by using Zonal polynomials techniques 
[32]. 

p(xi,x 2 ) = ^-e-^^u\{ Xl -x 2 f)e-^-*^\x 1 -x 2 \ (22) 
Ztt Z 

where I is the Bessel function of complex argument. 

In conclusion, we have expressed the joint distribution for the eigenvalues of our model 
in terms of Jack polynomials but we have not succeeded to derive explicit expressions for 
the correlation functions. 



3 Critical spectral kernel and density-density corre- 
lations of the C-S model at finite temperature 

We recall that the two level spectral function of our model is identical to the density- 
density correlation function of the C-S model (13) for A = 1. In [16] it was conjectured 
that the low temperature limit of the connected density-density correlation function of 
the C-S model at A = 1 is given by, 

( P (x)p(0)) T - (p(x)) T (p(0)) T = Rl c (x,0) = -K*(x,0) - Q^rfoO)) jH K T (t, 0)(23) 
where 

^> / ^ sin(7ra;) 

K T (x,0)=T , 1 > 24 
smn(7nri ) 

is the kernel of the C-S model (13) for A = 2. The temperature T = ^ and h is related 
to t and oo through (17) and (9). This result is valid in a normalization such that the 
average density of the particles is equal to unity. 

1 Some interesting results for for small values of N were obtained in [33]. 
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The idea is that, based on the Luttinger liquid nature of the C-S model [36], the known 
relation between the density- density correlation of the C-S for A = 1 at zero temperature 
and the spectral correlations of the GOE can be extended to finite low temperature. 
One simply replaces the kernel at zero temperature, which physically corresponds to free 
fermions for all invariant RMT ensembles 2 , by its finite temperature analogue [2] given 
by, 

^)=E^|g (25) 

where ip n are the single particle wave functions for free fermions and E n = ujn. The 
fugacity z for the free fermion distribution in (25) is determined by the total number of 
particles through, 

r 00 1 
N = / dxp(x) = Y . (26) 

where p{x) = K T (x,x) is the average spectral density. In the low temperature limit, 
ujt <C 1, we find 

* _1 = w 

In eq. (9) the quantities r and uo have been related to the parameter h of the matrix 
model (1). In the large N limit these relations simplify to 

" T ~m (28 > 

u ~ 2hN. 

and for h <C 1 the fugacity is given by z ~ e l l h . To obtain the average particle density for 
i^Owe can approximate the single particle wave function by plane waves with energy 
given by A; 2 . 



Pf) = if*r^np~^ (29) 



1 f°° „ 1 /K N "/h 

- / -i i n Tk 2 ~ V2 — 

IT JO 1 + Z 1 e TK 7T 

The unfolded spectral kernel is thus given by, 

= r(xM0),0) = ^ eo^Wg) 1 
V ; p(0) Jo 2v^ l + ^e* V ; 

In the low temperature limit h <C 1, the above spectral kernel coincides with (24) for 



irh 



T - 

± 2 



2 For A ^ 2 the particles obey exclusion statistics 
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The number variance of the eigenvalues near the center of the band is obtained by 
integrating the two point connected correlation function E% c (s, 0) including the self- 
correlations (23), 

£ 2 (L) = L + 2 ( L/m ds(L - s)R T 2 = 7Th/2 {s, 0). (31) 



2.5 




Figure 1: The number variance S 2 (L) (31) versus L for h = 0.1, h = 0.15 and h = 0.2. 
The linear behavior of the number variance is a signature of critical statistics. The slope 
for h <C 1 is x = h. 



The number variance S 2 (L) measures the stiffness of the spectrum. The fluctuations 
are small for the GOE with S 2 (L) proportional to log(L) for L > 1. For the Poisson 
ensemble, which is an ensemble of diagonal random matrices, the eigenvalues are un- 
corrected and £ 2 (£) = L. For critical statistics the number variance is asymptotically 
proportional to For % << 1 the slope has been connected with the multifractal dimen- 
sion D 2 observed in the wave functions of a disordered system undergoing a localization- 
delocalization transition [34, 35], 



where d is the spatial dimension of the system to be studied. As observed in Figure 1, 
the number variance of our model is linear for L ^> 1, with a slope x = h for h 1. 
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This linear behavior together with the absence of subleading logarithmic terms in the 
asymptotic behavior of the number variance suggests that our matrix model describes 
critical statistics. 



4 Testing the Kravtsov-Tsvelik Conjecture 

Below we discuss two independent methods to test the Kravtsov-Tsvelik conjecture. 
4.1 Conformal Calculation 

We review first [36, 16] how conformal techniques can be utilized to calculate the low 
temperature large distance asymptotic behavior of the two-point correlation function of 
the C-S model. 

Conformal field techniques [36, 37, 38, 39] can be used to compute the asymptotic 
behavior of the correlation functions of 1 + 1 dimensional systems with a linear gapless 
spectrum in the limit of large N number of particles, constant density n = N/ L and low 
temperature. 

In order to identify the conformal theory associated to the low energy excitations of the 
1 + 1 dimensional system one needs the value of the conformal anomaly c of the associated 
conformal field theory. Usually, c is obtained from the leading low temperature behavior 
of the free energy of the system. Then, the conformal weights of the primary fields of the 
conformal theory must match the leading low energy excitations of the 1 + 1 dimensional 
quantum system. The latter is usually evaluated either numerically or by finite size 
scaling and Bethe ansatz techniques. This program was carried out for the C-S model by 
Kawakami and Yang [36] . They found that the low energy collective excitations of the C-S 
model are described by two quantum numbers: AN related to excitations that change 
the number of particles and AD associated with excitations that move one particle from 
one Fermi point to the other. The energy, E, and momentum, P, of the leading finite size 
excitations of the C-S model are given by, 



E = AD + AN 2 + — i '-, 

LA LA L 

P = 2nAD + P L , 
2n 

P L = —[AN AD + N + + N~], (33) 

Lj 

where Pl stands for the momentum of the finite size excitations and iV + and iV~ label 
the conformal towers of states (secondary fields) associated to each primary field. Finally, 
they argued, based on thermodynamics arguments, that the conformal anomaly associated 
with the C-S model is c = 1. 
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In a system with conformal symmetry, the eigenvalues of both the Hamiltonian and 
the momentum are related to the right/left conformal weights x ntm /x n ^ m of the primary 
fields through the following relation, 

L , 

X n ,m = —[E + P L \, 

Xn,m = ^[E-P L ], (34) 

where we assume that E and Pl depend on quantum numbers n and m. Therefore, if 
we manage to find a conformal field theory with c = 1 and eigenvalues of the momentum 
operator and energy operator given by (33), the correlation functions of the C-S model in 
the asymptotic limit can be easily evaluated by means of conformal techniques. It turns 
out that the simplest conformal model with such properties is a free boson compactified 
on a circle with radius R = 1/ y/\. 

In general, observables do not have definite conformal dimensions and must be ex- 
pressed as a linear combination of conformal excitations. Since such conformal fields only 
describe the excitations close to the ground state one first has to decompose the expan- 
sion of observables into "fast" and "slow" modes [41]. The "slow" modes are described 
by the conformal fields and the the fast "ones" correspond to momenta that remain finite 
in the thermodynamic limit, i.e. to excitations with AD ^ 0. The density operator can 
be expanded as 

oo 

p(x) T = c m , n e i2 ^> n , , m (x) T (35) 

m,n=— oo 

where ip n ,o,m( x ,Q)T stands for the primary state (AD = n,AN = 0) associated with 
the above conformal theory, the phase in the expansion represents the momentum of the 
ground state for L — > oo (fast mode) and the index m accounts for the contribution of 
secondary fields. Since the density operator does not change the number of particles only 
excitations with AiV = contribute to the expansion. The coefficients c m ^ n are found 
from the zero temperature limit (GOE). 

The density-density correlations at finite temperature can be easily obtained using the 
known result for the correlation functions of the conformal fields ^> n ,o,m(x,0) T - The first 
terms of the conformal prediction for the density-density correlations of the C-S model at 
A = 1 in the low temperature limit are thus given by 

Rl(x, 0) = (p(x)p(0)) T ^ - — cos [ 2nx \ _ 3 ^ (36) 

sinh (ttxT) 2 sinh (irxT) 2 sinh (ttTx) 

With a temperature as given by the Kravtsov-Tsvelik conjecture, i.e. T — the confor- 
mal result coincides with the asymptotic expansion of the conjectured result (23). Since 
both results have been obtained by using completely different methods, this calculation 
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supports the validity of the kernel (24) in the low temperature, long distance limit. 
4.2 Susceptibility 

The slope of the large distance asymptotic behavior of the number variance is determined 
by the isothermal susceptibility of the C-S model which can be obtained from the C-S 
partition function. On the other hand, this slope is determined by an integral of the 
unfolded two-point cluster function Y 2 (r) according to 

/oo 
Y 2 (r)dr). (37) 
-oo 

Therefore, agreement between the conformal calculation and the Kratsov-Tsvelik conjec- 
ture for the large distance asymptotic behavior of the two-point correlation function does 
not necessarily imply that the asymptotic behavior of the number variance is given by 
the susceptibility. 

The susceptibility x m t ne grand canonical ensemble which measures the fluctuations 
of the number of particles in a box of length L, 

X = (N 2 ) - (AT) 2 , (38) 

can be expressed as 

X = ^(N), (39) 

where Z\ stands for the fugacity and (N) is the average number of particles. Remarkably, 
the C-S gas for arbitrary statistical coupling A can still be considered a free gas but 
with exclusion statistics [31]. Indeed, Sutherland [19] has shown, by using the Bethe 
Ansatz and a method previously developed by Yang and Yang [40] for the Bose gas with 
a delta interaction, that the occupation number n(k) of a gas of C-S particles satisfies the 
following transcendental equation, 

(1 - An(A;)/2) A / 2 (l + (1 - \/2)n{k)f- x / 2 = n(k)e Te ^/z v (40) 

For the special case A = 1, corresponding to the GOE, we obtain 

n(k) = \ (41) 

1 ' ^1 + ^e 2 -« 1 ' 

where z\ is the fugacity for this distribution function and e(k) = k 2 is the energy of a 
single particle. To find the relation between the fugacity and the parameter h in the 
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Moshe-Neuberger-Shapiro model for A = 1 we have to use the single particle energies 
corresponding to the C-S model (12). We thus have 

1 roo roo 1 

N = — dx dk , (42) 

7T J-oo J-oo U _|_ ^ e 2T(k 2 + \u) 2 x 2 ) 



The asymptotic behavior for large z 1 can be obtained easily by changing to polar coordi- 
nates. This results in 

LUT 

Using (9) we then find in the limit of small h, 

Zl = e^ 2h . (44) 

Since the density of particles is x-dependent in a harmonic box we calculate the sus- 
ceptibility for particles in a rectangular box with fugacity given by (44). In this way the 
susceptibility can be compared to the slope of the number variance which is calculated in 
the center of the spectrum. 

The average number of particles in a box of length L is given by 

w "L d + vr (45) 



If the fugacity is parameterized as z'f = e 2Tfe2 we have in the low-temperature limit, 

(N) = 2k-. (46) 

7T 



Then 



dzi 71 Jo (1 + 4^ 



After the change of variable 5k = k — k and expanding around the Fermi surface we find, 



L f°° 
X = - d5k 

7T J-oo 



8e 



Ark&k 



oo (1 + 4 e 4Tfc5fc^ 

(N) (N) 



2rk 2 21og£i 



= h(N). (47) 



The above result should be compared with the calculation of the asymptotic behavior 
of the number variance from the two-point spectral correlation function, 

X ~£ 2 ((iV» for (JV)-oo, (48) 
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where S 2 ((iV)) is denned by 

r {N) 

S 2 ((iV)) = (N 2 ) - (iV> 2 = (N)-2 dr((N) - r)Y 2 (r). (49) 

Jo 

Here, Y 2 {r) is the unfolded two-point cluster function. According to the Kravtsov-Tsvelik 
conjecture it is given by 

Y 2 {r) = K 2 (r) + f°° K(r')dr', (50) 

where K{r) is the kernel 



oo 



cos(Ax/p) 



K(r) = / dk V^j, (51) 



and p is the average spectral density. For this cluster function, and, in fact any cluster 
function that decreases stronger than 1/r, we recover the relation (37) for (N) — > oo. 
With the fugacity parameterized by z 2 = e rfc2 = e l ^ h (see section 3) we have in the 
low-temperature limit 

K(0) = 1 = — . (52) 
np 

After partial integration of the second term of (50) and using that K(0) = 1 we obtain, 

E 2 ((N)) = (N) -2(N) / (2K 2 (r) - K(r))dr + 0(N°). (53) 

Jo 

Integrating by parts and making an expansion about the Fermi surface k results in, 



sinkr/p r°° cos sr/2rkp 1 sin(kr/p) 
4nr J-oo cosh 2 s/2 2kpr sinh(7rr/2A;pr) ' 

The integral over r in (54) can now be performed analytically resulting in the suscep- 
tibility 



x= w = <m = h(Ny (55) 

r k 2 log z 2 

This slope is in agreement with the result obtained from the partition function of the C-S 
model. In agreement with our naive expectation, the value of the slope is a factor 2 larger 
than the one found for the original Moshe-Neuberger-Shapiro model for A = 2 [2]. 

In the kernel (51) the momentum integral is weighted by the occupation number which 
in this case is the Fermi-Dirac distribution. Since the occupation number of the C-S model 
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for A = 1 is given by (41) it seems more natural to make this choice instead. This results 
in the kernel 



where z 1 = e^ T = e 1 ^ 1 ' is the fugacity. We choose the normalization of the kernel such 
that K(0) = 1. Then the zero temperature limit (r — > oo) of this kernel is the usual 
sine-kernel, smnx/nx. Below we show that the conjecture (56) disagrees with both the 
conformal calculation and the susceptibility (47). 

Let us first derive the large x asymptotic behavior of the kernel (56). In the low 
temperature limit the average number of particles is again given by (N) = 2Lk/irp, with 
normalization condition 2k/irp = 1. After partial integration the integral can be rewritten 
as 

_ 1_ roc grkz^ 2 sMZkx/p) 
{ ) Tix Jo (l + 4^ 2 e 2 ^ 2 )3/2 • 

In the low temperature limit the integrand is strongly peaked at k ~ k, and the integral 
can be calculated by a steepest descent approximation 

2 I" 00 /^ e 2ix(k+u/Tk)/p e 4u 

K(x) ~ Im — / du- 

71X J-, 



1 +4e 4 «) 3 /2 



I m f 71 ^ \ e 2ikx/p-nxk/2T-nix(\og4)/2Tk ^g-j 



37T V rk J 



This asymptotic result is in disagreement with the prediction from the conformal calcu- 
lation. 

Next we compare the asymptotic behavior of the number variance with the suscepti- 
bility. The comparison of the susceptibility can again be made by computing the asymp- 
totical behavior of the number variance. For an exponentially decreasing kernel we have 
previously shown that the asymptotic behavior of the number variance is given by 

E 2 ((AO) = (N) - 2(N) / (2K\x) - K{x))dx + 0((N) ). (59) 

Jo 

where K(x) is the kernel (56) with average eigenvalue spacing normalized to unity. Using 
that 

f°° sin irax 1 

/ dx = -, 

Jo nx 2 

oo g4s I 

00^(1 + 4^)3/2 = g. («) 
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one easily shows that in the low temperature limit, 

roo i 

J o K(x)dx = -. (61) 
The other integral in (59) can be written as 

roo roo 128 fV -11 - 1l' 

/ K\x)dx= dx——^ dy dy'sm((k + ^ r )2x/p)sm((k + ^ T )2x/p)G(y)G(y'). 

JO JO {TTX) Z J-oo J-oo tk rk 



(62) 



where 

-Ay 



e 



The integral over x can be evaluated using the formula 

r °° sin (ax) sin (ox) an 



This results in 



f°° , sin (ax) sin (ox) an . . 

/ dx — ^ — \ v ; = — for a < 6. 64 

Jo x 2 2 y J 

rOO j roo ry rji 1 

/ K\r)dr = - + 128 / dy dy>^L-_G(yW). (65) 
Jo 2 J-oo J-oo nrkp 



nrkp 

The asymptotic behavior of the number variance is thus given by 



= -64 



(N)_ 
j 

/CO 
. -co dy (T+ 4e 4 f ) 2 

2(N) 



£ 2 ((iV» = -512^ J°°Jy f jy'y'G{y)G{yi) 
(N) [°° ye Ay 



nrkp 

(N) 
log Zj 



log 2 

log 2. (66) 



To obtain the expression after the second equality sign we have performed a partial 
integration using the identity 

G(y) = --4- , 1 :■ (67) 

yU> 8 dy VI + 4e 4 ^ 1 ; 

This result for the susceptibility differs from the result obtained from the thermodynamic 
properties of the C-S gas. We conclude that the kernel (56) does not describe the corre- 
lations of the critical Random Matrix Model (1). 
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5 Critical statistics and quantum chaos 



In this section we introduce the concept of multifractal wave functions in the context of 
the Anderson transition and show how it may be relevant in the study of deterministic 
quantum chaotic systems. 

By now it has been well established that the appearance of critical statistics at the 
Anderson transition is intimately related with the multifractal properties of the wave- 
functions [9, 34, 35, 43]. We wish first to introduce intuitively the concept of multifractal 
wave functions [42] . 

Let us consider the volume of the subset of a box for which the absolute value of the 
wave function \1/ is larger than a fixed number M. If this volume scales as L d * (with 
d* < d), then d* is called the fractal dimension d* < d of In case the fractal dimension 
depends on the value of M, the wave function is said to be multifractal. More formally, 
multifractality is defined through the inverse participation ratio, 



where \l/ n is the wave function with energy E n , and D p < d is a set of exponents char- 
acterizing the anomalous (multifractal) scaling of the moments of the wave function. We 
remark that, although confined to fractal subsets of the sample, wave functions of such 
systems overlap strongly when their energies are close enough [43]. Such strong overlap 
is responsible for the short-range level repulsion observed at the Anderson transition. 
It is worthwhile to note that this anomalous scaling has, in principle, a pure quantum 
mechanical origin. As the density of impurities increases, the deBroglie wavelength of 
the particles becomes comparable with the mean free path and localization effects start 
to be relevant. We stress that the classical dynamics of the Anderson transition does 
not provide us with valuable information to describe quantum spectral correlations. One 
may wonder to what extent such multifractal behavior may be observed in deterministic 
quantum chaotic systems. 

What has become known as the Bohigas-Giannoni-Schmit conjecture [44] is that 
generically quantum spectra of classically chaotic systems are correlated according to 
the Wigner-Dyson random matrix ensembles, whereas spectral correlations of classically 
integrable systems are close to Poisson statistics. In most cases, by modifying the param- 
eters of the system, a transition from integrable to chaotic dynamics can be observed. If 
the KAM theorem is applicable, this transition is smooth and both integrable and chaotic 
regions coexist until the last KAM torus is completely destroyed. Although spectral 
statistics of such mixed system have been described in terms of banded random matrix 
models [45], they are believed to be non-generic and different form critical statistics [46]. 




(68) 



r 
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Figure 2: The spectral rigidity A3 of the anisotropic Kepler problem obtained in [52] 
versus the prediction of our model for h = 0.16 obtained from (30). Both curves are 
barely distinguishable. The numerical data are reprinted from Fig 2. in Ref. [52] 



The situation is different in cases where the KAM theorem does not apply. In those 
systems, the invariant KAM curves may not exist at all and small changes in the cou- 
pling constant can produce qualitative modifications in the classical phase space. The 
dynamics is, in general, intermediate between chaotic and integrable. The lack of KAM 
tori permits a particle explore the full available classical phase space without having full 
chaotic motion. Such forms of phase space are also known as stochastic webs [47]. In 
certain cases, the classical phase space becomes increasingly intricate, showing both self- 
similar and fractal properties [47]. We notice that such a structure is reminiscent of the 
way that the KAM tori break up into "fractal" orbits of zero dimension [48] (cantori) 
as the system becomes chaotic. Classically, cantori represent strong obstacles to phase 
space transport. Our aim is to study the effect, if any, of such self-similar structure in 
the spectral correlations of the quantum counterpart. Roughly speaking, the influence 
of cantori on the quantum dynamics will depend on the relation between the size of the 
cantori and Planck's constant. For cantori smaller than the Planck cell, quantum dynam- 
ics cannot resolve the classical fine structure. In this case, cantori act as perfect barriers 
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to the quantum motion resembling the effect of a classically integrable system and the 
spectral correlations of the quantum counterpart are close to Poisson statistics. In the 
intermediate case the situation is less clear. Recently, it has been reported [49] that can- 
tori drive spectral correlations smoothly from Poisson to RMT as the system approaches 
to the ergodic regime. 

Below, we present numerical evidence that deviations from GOE statistics caused by 
the self similar structure of the classical phase space may be described by critical statistics 
at least while the deviations from GOE are small. 

5.1 The Anisotropic Kepler problem 

The anisotropic Kepler Hamiltonian 

H = \p 2 p + \ivl - \ (69) 

is an interesting example of a non-KAM system undergoing an abrupt chaotic integrable 
transition. It has been utilized as a model of donor impurities in a semiconductor [50, 51]. 
Even for small departures from the integrable case, 7 = 1, the classical phase space is 
densely filled with remnants of cantori [52]. Gutzwiler has shown that for 7 < 8/9 the 
orbits in phase space can be uniquely represented in terms of symbolic dynamics. Such 
representation is a signature of hard chaos. Indeed, for 7 < 1/2 there are no islands 
of stability in phase space. Furthermore, the measure of the surface of section based 
on the symbolic dynamics is multifractal with respect to the usual Liouville measure. 
Since the periodic orbits can be effectively enumerated, the energy levels of the quantum 
counterpart can be approximately evaluated by means of analytical techniques [53]. 

A numerical study of the spectral correlations of highly excited states of this system 
was carried out in [52]. In a basis in which the Hamiltonian has a band structure, they 
succeeded to obtain up to 5500 energy levels. In Fig. 2 we show their result for the 
spectral rigidity of the spectrum from level 2501 to 5500. As observed, the deviations 
from the GOE are very well described by the critical random matrix model (1). Based on 
the analogy with disordered systems, we conjecture that the wavefunctions of this system 
are multifractal. We are not aware of numerical results that can confirm or disprove this 
conjecture. 

5.2 Kicked particle in a infinite potential well 

Recently, in [54], another non-KAM system with similar properties, a kicked particle in 
a infinite potential well, was studied both quantum mechanically and classically. The 



20 




20 L 40 60 



Figure 3: The number variance £ 2 (£) versus L. The number variance of the quasi-energy 
levels of a kicked particle in an infinite potential well at K = kT = 50 [54] agrees with 
the prediction of the critical GOE (31) for h = 0.2122 (upper curve) up to 10 eigenvalues. 
The downward tendency of the numerical result may be due to finite size effects. The 
error in the numerical results is indicated by the thickness of the curve. We thank to 
Baowen Li [54] for kindly providing us with 1024 energy levels to compute the number 
variance. 

Hamiltonian is given by, 

H = y — + V(x) + kcos( y x + l) S{t-nT) (70) 

n=— oo 

where V(x) is an infinite well potential of length ir, T is the period of the kick and k the 
strength. Concerning the classical motion, the KAM theorem is not applicable because 
the potential is not smooth. Indeed, it was found [54] that the classical phase space 
resembles a stochastic web with a self similar structure. This is in contrast with the 
standard kicked rotor where the classical phase space is a mixture of chaotic an integrable 
parts separated by KAM tori. 

The quantum mechanical properties of the model are described by the evolution op- 
erator U = e - t P 2T / A e kcos ( x + l ) e - l P 2T l i over a period T of the kick. The quasi energies 
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associated to this operator were obtained in [54] by diagonalizing U in a basis of 1024 
eigenstates of the free Hamiltonian (q\n) = sin(ng). 

Unlike the kicked rotor where the matrix evolution has an exponential decay in a basis 
of plane waves, it can be shown that the matrix elements of U are well described by a 
random banded matrix (RBM) with power-like decay, \{m\U\m + n) \ oc ^ for b <C n and 
constant for b ^> n where b ~ K = kT is the size of the band. 

In agreement with the results of [3], the nearest neighbor distribution reported in 
[54] smoothly interpolates between Poisson and GOE as K is increased. Recently, an 
experimental realization of this model was studied in [55]. 




Figure 4: The spectral rigidity A 3 (L) obtained in [54] for the energy levels of a kicked 
particle in an infinite potential well at K = kT = 50 is accurately described by the critical 
kernel (30) at h = 0.2122 (upper curve). The error in the numerical results is indicated 
by the thickness of the curve. 

Analysis of results 

In Figures 3 and 4, we show the number variance and the A 3 -statistic of the sequence of 
1024 eigenvalues obtained in [54] for kT = 50. The upper curve is the analytical result 
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derived from the two-point correlation function (23) with kernel (30). In both cases the 
numerical result is plotted with its error (see below). 

For the A 3 -statistic, the x 2 is minimized on the interval [0, 30] for h = 0.2122 with a 
value 3 of x 2 — 0-32. The errors in the definition of x 2 h &ve been calculated by splitting 
the 1024 eigenvalues into eight ensembles of 128 eigenvalues and evaluating the number 
variance for each ensemble separately (denoted by E 2 (L), % = 1, • • • , 8). The error in the 
number variance is thus given by 

where S mcan (L) = i ]T £ 2 (L). (71) 
8 i=i 

The error in the A3-statistic is obtained from this error by means of a Monte-Carlo 
simulation using the relation between A 3 (L) and S 2 (L) [56]. As observed in Figure 4, the 
error in A 3 (L) is much smaller that the error in E 2 (L). 

Next we ask the question whether the asymptotic behavior of the spectral rigidity is 
linear without a logarithmic correction, just as in the analytical case. If a logarithmic 
correction is absent, one can almost discard the possibility of a mixed classical phase space 
as the reason of the observed deviation from the GOE. In order to prove the absence of 
such term we fit the numerical curve A 3 (L) num to A 3 (L) fit = a + bL + clogL. For instance, 
on a interval [15,25] a best fit is obtained for a = 0.110 ± 0.004 , b = 0.0156 ± 0.0001, 
c = 0.002 ± 0.002 with a value of x 2 = 0.015. We find that the value of the coefficient c 
is compatible with zero. This suggests that the classical phase space is not a mixture of 
chaotic and integrable regions. 

In the case of the number variance, because of the size of the error, no conclusive 
evidence on the absence of the logarithmic term can be obtained by such fit. At large 
distances the number variance seems to deviate from a linear behavior by a quadratic 
term. Although such terms are typically caused by finite size effects, we do not have a clear 
understanding of its origin. Since the A 3 -statistic projects out a quadratic dependence 
of the number variance the linear behavior persists to much larger distances in this case 
(see Figure 4). 

Finally, let us confront the conjecture (23) with our present numerical results. As we 
mentioned previously, it may look reasonable to replace the kernel (30) by (56). It can be 
shown that the spectral rigidity obtained from the kernel (56) is almost indistinguishable 
from the one obtained from (30) and therefore the agreement with the numerical result is 
expected to be equally good with one fitting parameter at our disposal. However, a more 
careful analysis shows that the kernel (56) leads to a value of x 2 much higher than the 
one obtained from (30). For the interval [0, 30], using the kernels (30) and (56), a best fit 

3 Here and below we find values of x 2 that are significantly less than one. This is possible because the 
values of A 3 (L) for different values of L are correlated. Therefore, our values of \ 2 have to be used with 
care and cannot be interpreted in terms of a x 2 distribution. 



a(L) 



V8 



i=l 
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is obtained for h = 0.2122 with a value of xlermi = 0-32 and for h' = 0.325 with a value 
of Xlqrt = 2 - 27 ? respectively. 

The above findings suggest that a self-similar classical phase space dominated by 
cantori has a strong impact on the quantum spectral correlations. Critical statistics 
appears as the leading candidate to describe such correlations and, consequently, enlarge 
the range of applicability of random matrix ensembles. 

Finally, we list other quantum systems between integrable and chaotic whose quantum 
spectral correlations show a similarity with critical statistics: quantum billiards with a 
point scatterer [57, 58], the Kepler billiard [17, 59], semiconductor billiards [61, 62], the 
stadium billiard inside certain range of parameters [63]. For applications concerning 
pseudo integrable billiards we refer to [17]. 

6 Conclusions 

In this article we have introduced a one parameter ensemble of symmetric random ma- 
trices. This ensemble interpolates between the Gaussian Orthogonal Ensemble and the 
Poisson ensemble and is capable of describing critical statistics. 

We have shown that, in an eigenvalue basis, the joint eigenvalue distribution of our 
model coincides with the diagonal density matrix of the C-S model at finite temperature 
where the additional parameter of the matrix model plays the role of temperature. Re- 
markably, this equivalence can be extended to all random matrix ensembles associated 
with the large families of symmetric spaces according to the Cartan classification thus 
providing a novel link between strongly interacting quantum systems and random matrix 
theory. 

We have calculated the spectral correlation functions based on a recent conjecture by 
Kravtsov and Tsvelik for the correlation functions of the C-S model in the low temperature 
limit. We have tested the validity of this conjecture by two independent methods: one 
based on the effective conformal symmetry of the C-S model in the low temperature limit 
and the other based on the thermodynamical properties of a gas of particles governed 
by the C-S Hamiltonian. We have found that both the long distance low temperature 
behavior of the two-point correlation function obtained from the conformal calculation 
and the susceptibility of the C-S model agree with the conjecture made by Kravtsov and 
Tsvelik. 

Based on the Kravtsov-Tsvelik conjecture we find that, although level repulsion is 
still present, the number variance is asymptotically linear with a slope less than one and 
no sub leading logarithmic term present. This indicates that our random matrix model 
describes critical statistics. 

Finally, we have argued that critical statistics is relevant to describe spectral corre- 
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lations of chaotic quantum systems for which the Poincare section of the classical coun- 
terpart is globally self-similar/fractal. Two examples with such classical phase space, a 
kicked particle in a potential well and the anisotropic Kepler problem, have been discussed 
in detail. In both cases, long range spectral correlators such as the number variance and 
the A 3 -statistic are accurately described by our analytical results based on the Kravtsov- 
Tsvelik conjecture. Indeed, for the kicked particle, we have shown that the spectral 
rigidity is asymptotically linear with no subleading logarithmic term present. This may 
be an indication that the wavefunctions of this model show multifractal properties. 
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